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ABSTRACT 

The Initial Mass Function (IMF) of early-type galaxies (ETGs) has been found to 
feature systematic variations by both dynamical and spectroscopic studies. In par¬ 
ticular, spectral line strengths, based on gravity-sensitive features, suggest an excess 
of low-mass stars in massive ETGs, i.e. a bottom-heavy IMF. The physical drivers 
of IMF variations are currently unknown. The abundance ratio of a elements, such 
as [Mg/Fe], has been suggested as a possible driver of the IMF changes, although 
dynamical constraints do not support this claim. In this letter, we take advantage of 
the large SDSS database. Our sample comprises 24, 781 high-quality spectra, cover¬ 
ing a large range in velocity dispersion (100 <<jq< 320 kms -1 ) and abundance ratio 
(—0.1 <[Mg/Fe]< +0.4). The large volume of data allows us to stack the spectra at 
fixed values of ao and [Mg/Fe]. Our analysis - based on gravity-sensitive line strengths 
- gives a strong correlation with central velocity dispersion and a negligible variation 
with [Mg/Fe] at fixed ao- This result is robust against individual elemental abundance 
variations, and seems not to raise any apparent inconsistency with the alternative 
method based on galaxy dynamics. 

Key words: galaxies: stellar content - galaxies: fundamental parameters - galaxies: 
formation 


1 INTRODUCTION 


The stellar IMF, i.e. the mass distribution of stars in a stellar 
population (hereafter SP) at birth, has been largely assumed 
to be universal, mostly because of lack of evidence regarding 
variations of the IMF a mong star clusters and OB associa - 
tions in our Galaxy (e. g. I C h ab r i e rll+Xfl : I K r o up a et al.ll2010l) . 
However, recent studie s, based on independent techniques, 


201a 

: ITortora et al. 


20ia 

), lensing (e.g. Auger et al. 

201c 

; but see Smith & Lucev 20131), and spectroscopy 

(e.g. [Cenarro et al. 1 

20031; van Dokkum & Conrov 

201C; 


Spiniello et al. 


La Barbera et al. 


2014; 

Horn 


as well as F 
hereafter F13 


2013 


and LB13, respec¬ 
tively), have found that the IMF varies systematically with 
galaxy mass in ETGs. 

Th e origin of t his systematic var iation remains un¬ 
known. 1 Conroy fc van Dokkuml d2012bl . hereafter CvD12b) 
found, for a sample of 38 nearby ETGs, that the IMF nor¬ 
malization (i.e. the stellar M/L, normalized with respect to 


a MW-like IMF) correlates more strongly with the abun¬ 
dance of a-elements, such as [Mg/Fe], than with central 
velocity dispersion (ao), claiming that [Mg/Fel, w hic h is a 
proxy for the star-form ation time-scale dThomas et al.ll2005l : 
Ide la Rosa et al]l201ll ). might be the main driver of the IMF 
in ETGs. However, this finding is challenged by the fact 
that dynamical constraints, for a subsample of 34 ETGs in 
comm on between CvD12b and ATLAS 30 llCap pellarietall 
I2013I) . correlate better with ao than [Mg/Fe] ( Smithl 2014 ). 
Furthermore, no correlation of dynamical constraints with 
SP parameters has been d etected in the ATLAS 3D sam¬ 
ple (IMcDe rmid et alj 120141 ). Given the importance of con- 
structi ng a com prehensiv e theory o f a non-universal IMF 
(e. g. [Hopkin sj 120131 : Ichabrier, Hennebelle, Chariot! 12014 
iFerreras et al.ll2015t ). the above inconsistencies call for fur¬ 
ther investigations, based on larger samples, of the possible 
drivers of IMF variations. In LB13, we assembled a battery 
of 126 stacked spectra, with varying ao and [Mg/Fe], finding 
that gravity-sensitive features in ETGs depend very mildly 
on [Mg/Fe] at fixed a 0 . Such a result suggests no correlation 
between IMF and [Mg/Fe]. In the present letter, we comple- 
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ment the results of LB13, performing a quantitative analysis 
of how the IMF slope depends on both <jq and [Mg/Fe]. 


2 DATA 

The SPIDER sample llLa Barbera et ahl l2010al . hereafter 
LB10) consists of 39,993 nearby (0.05 < z < 0.095) 
ETGs, selected fro m Data Release 6 of the Sloan Digital 
Sky Survey ( SPSS. l~Adelman-McCarthy et alJl2008f ). ETGs 
are defined as bulge-dominated systems, featuring pas¬ 
sive spectra within the SDSS fibres. All SPIDER. ETGs 
have SDSS spectra available, covering the spectral range 
from 3800 to 9200 A, with a spectral resolution Ui n3 t ~ 
60 kms -1 . For each spectrum, we use its central veloc¬ 
ity dispersion, ao, from the SDSS database, and an em¬ 
pirical proxy, [ZMg/ZFe], for the [Mg/Fe] abundance ra¬ 
tio (see LB13 for details). The [ZMg/ZFe] is the differ¬ 
ence of total metallicity estimates, [Z/H], measured from 
the Mgb5177 and Fe3= (Fe4383 + Fe5270 + Fe5335)/3 
spectral indices, respectively. Total metallicities are mea¬ 
sured at fixed age (estimated with the spect ral fitting 
code STARLIGHT; ICid Fernandes et alj 1200511 , compar¬ 
ing observed line strengths to predi ctions of MILE S Sim - 
ple Stellar Population (SSP) models dVazdekis et al.lf2012lb 
As shown in LB13, [ZMg/ZFe] is tightly correlated with 
[a/Fe], derive d from SP models taking abun dance ratios 
into account (IThomas. Maraston fe Johanssonl 120111 1. with 
[a/Fe] ~ 0.55 [ZMg/ZFe] - Since Mg is the reference a-element 
of an SP model, and both Mgb5177 and Fe3 are mostly 
sensi tive to the abundance of Mg and Fe, respectively (see, 
e.g.. [Johansson et al.|]201^ . hereafter JTM12), in this let¬ 
ter we use [Mg/Fe] rather than [a/Fe], with the conversion 
[Mg/Fe]~ 0.55[ZM g /ZFe]. Notice that since [ZMg/ZFe] actu¬ 
ally measures [Mg/Fe], our conclusions do not rely on the as- 
sumption that a elemen ts should track each other (see, e.g., 
iKuntschner et ad] l2010l l. Following LB13 and F13, we re¬ 
strict the analysis to a subsample of 24,781 SPIDER ETGs, 
with (i) 100 ^ (Jo ^ 320 kms -1 , (ii) low internal extinction 
(i.e. a color excess estimate E(B — V) < 0.1), and (iii) better 
signal-to-noise (S/N) ratio (see below). We classify the spec¬ 
tra into 18 narrow bins of ao(see LB13 for details), each bin 
with a width of 10 kms -1 , except for the last two, for which 
we adopt the range [260, 280], and [280, 320] kms -1 , respec¬ 
tively. For each bin, we only select spectra whose S/N ratio 
is higher than the lowest quartile of the S/N distribution in 
that bin. 

In the present letter, we analyze the median-stacked 
spectra assembled by LB13 for this sample, i.e. (i) the 18 
“(Jo stacks”, obtained by selecting all ETGs in each ao bin, 
and (ii) the “[ZMg/ZFe] sub-stacks”, assembled by selecting, 
for each bin in ao, those ETGs (1) below the 10-th, (2) below 
the 25-th, (3) between the 25-th and 50-th, (4) between the 
50-th and 75-th, (5) above the 75-th, and (6) above the 90- 
th percentiles of the [ZMg/ZFe] distribution. This procedure 
yields 18 ao stacks, plus 18 x 6 = 108 [ZMg/ZFe] sub-stacks, 
for a total of 126 stacked spectra. At fixed ao, the [ZMg/ZFe] 
stacks span a range of d([ZMg/ZFe]) ~ 0.6 (corresponding to 
5([Mg/Fe]) ~ 0.3), allowing us to single out the effect of ao 
and abundance ratio on IMF slope. 


3 ANALYSIS 

We rely on the extende d MILES (MIUSCAT) SP mod¬ 
els (IVazdekis et alJl2012f l with age (total metallicity) 4 < t < 
14 Gy r (—0.4 <[Z/H[ < +0.22), and ten choices of a bimodal 
IMF dVazdekis et alj fl99fih ■ i.e. a power-law, with slope Ti, 
in log-mass units, smoothly tapered to a constant value to¬ 
wards low-mass stars (< 0.6M©), with 0.3 ^ Tj, ^ 3.3. 
Note that Pb~ 1.3 gives a good representation of a Kroupa 
(MW-like) distribution, while higher values of Pb correspond 
to bottom-heavier distributions. For each stacked spectrum, 
we derive Tb b y a similar proce dure to LB13 (see also F13 
and lMartm-Navarro et al.ll2015l . hereafter MN15). We min¬ 
imize the expression: 

X 2 (t,{ Z/H],P b ,[X/Fe]) = 

where Ei and Em,,, are observed and model line-strengths 
for a selected set of spectral features, [A/Fe] is the abun¬ 
dance ratio of different chemical elements (see below), Ai,x 
is the sensitivity of the i-th line strength to a given ele¬ 
mental abundance, i.e. Aj i x= <5(^M,i)/5([A'/Fe]), and ai 
is the uncertainty on Ei, obtained by adding in quadra¬ 
ture the statistical error on Ei with its intrinsic uncer¬ 
tainty (see sec. 6 of LB13). The Em,i are computed with 
MIUSC AT models, w hile the A,, x facto rs use the publicly 
available IConrov fe van Dokkuml (12012al . hereafter CvD12a) 
stellar population models, having solar metallicity, old age 
(t= 13.5 Gyr), and a Kroupa IMF. The free fitting parame¬ 
ters in Eq. |T| are age, metallicity, IMF slope, and elemental 
abundances. Both Em,i and A;_x are computed by matching 
the spectral resolution of the models to that of each stacked 
spectrum. Errors on best-fitting parameters are computed 
with a bootstrap approach, shifting observed line-strengths 
according to their uncertainties. Notice that by varying the 
Pb, we actually alter the relative fraction of low- to high- 
mass stars in the IMF by changing the behaviour at higher 
masses (> 0.6 Mg). Although this is different from changes 
in aunimodal (single power-law) distribution (e.g. CvD12a), 
in practice, gravity-sensitive features are only sensitive to 
the mass fraction of low-mass stars in the IMF (see figure 19 
of LB13). 

To test how abundance ratios may affect the estimate 
of Pb, we adopt two alternative methods, differing for the 
selected set of spectral features and elemental abundances 
included in the fitting (see Tab. [lj: 

1SSP (Single Burst) method- We use the same set of 
spectral features as in LB13, except for the Calcium triplet, 
CaT, which cannot be measured over the SDSS spectral 
range for most SPIDER ETGs, and is replaced by the Ca2 
feature alone (see MN15). No [X/Fe] term is included in the 
fits. Indeed, we found that, overall, the leading abundance 
ratios of gravity-sensitive features have a negligible impact 
on relative variations of Pb (see, e.g., green and cyan curves 
in figure 12 of LB13). In contrast to our previous works, 
we do not apply any [Mg/Fe] empirical correction to bring 
observed line-strengths to solar-scale, as (i) the variation 
of most gravity sensitive-features with ao is much stronger 
than that with [Mg/Fe] (LB13), and (ii) our aim, in the 
present letter, is to constrain the IMF as a function of both 
ao and [Mg/Fe], hence we do not want to apply, a priory, 
any [Mg/Fe] correction. 
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Figure 1. Trends of with ctq (left-) and [Mg/Fe] (right-panels) for our 126 SDSS-based stacked spectra. Upper and lower panels 
correspond to the 1SSP and lSSP+[X/Fe] fitting methods to infer the IMF slope (see Sec. 0. The <jo stacks (see Sec. 0 are plotted with 
filled circles, and connected by a black curve, while the sub-stacks with varying [Mg/Fe] are plotted with empty circles. Error bars are 
quoted at the 1 a confidence level. Blue-through-red colors (right) encode, as in LB13, the increase of cro; whereas green-through-orange 
(left) correspond to the range in [Mg/Fe]. Stacks with a large uncertainty (> 50%) on F|, are plotted in grey (bottom-panels; see the 
text). Notice the tight (poor) correlation of F|, with <tq ([Mg/Fe]) for both fitting methods. 


1SSP+[X/ Fe] method - We fit all individual abundance 
ratios, whose Ai,x can be computed with CvD12a mod¬ 
els, and adopt a much wider set of Lick-based features than 
in the 1SSP case, making sure to include, for each [X/Fe], 
at least one index with prominent sensitivity to it (based 
on figure 1 of JTM12), as summarized in Tab. [T] Notice 
that the values of A^x may not be constant as a function 
of age (as assumed in Eq. [T]) , as index responses to abun¬ 
dance variations could become milder toward s younger ages 
with respect to those found in massive ETGs llSansom et alj 
120131 ). implying that 1SSP models, instead of lSSP+[X/Fe], 
might describe better our 1ow-<t stacks (leaving unchanged 
the conclusions, see below). Moreover, since CvD12a mod¬ 
els (i.e. the Ai,x factors) work at fixed Fe, rather than total 
metallicity (MILES), the [Z/H] from Eq. |T] does not nec¬ 
essarily measure the total metallicity. This is nevertheless 
unimportant for the present letter, as we study trends with 
[ZMg/ZFe], rather than [Z/H] (see Martfn-Navarro et al., in 
preparation) . 

We notice that for each stack, the H/3 0 age indicator is cor¬ 
rected for nebular emission, as in LB13. 


4 RESULTS 

Fig- CD shows the best-fitting Tb, for all stacked spectra, vs. 
(Jo (left-) and [Mg/Fe] (right-panels), respectively. The 00 
([Mg/Fe]) stacks (see Sec[2]) are plotted with filled (empty) 
circles. Regardless of the fitting method (upper vs. lower 


panels), and in particular even when taking abundance ra¬ 
tios into account (lSSP+[X/Fe] method), we find a tight 
correlation between Tb and co, consistent with previous 
spectroscopic studies (e.g. CvD12b, LB13). Notice that the 
quality of the fits is generally good for all stacks, with a 
reduced x 2 ~ 2 (~ 1) for the 1SSP (lSSP+[X/Fe]) meth¬ 
ods (see sec. 7 of LB13), and - most importantly for the 
present work - no significant x 2 dependence on either ao 
or [Mg/Fe] is found. The Pb-(Jo relation is shallower for the 
lSSP+[X/Fe] case (see below), consistent with the fact that 
different Pb-uo slopes are found when adopting different sets 
of spectral features (F13. l^iniello_et_aL 20lJ][ i and/or dif¬ 
ferent methodologies (LB13. ISpiniello et al.ll2015lb For the 
present work, we do not aim to establish which method pro¬ 
vides a more reliable Tb-cJo relation, but, rather, to test the 
robustness of our results against two rather extreme cases, 
namely, the 1SSP method, where the impact of abundance 
ratios is neglected; and the lSSP+[X/Fe] approach, where 
we attempt to model the effect of individual elemental abun¬ 
dances, with theoretical (CvD12a) stellar population models 
(whose ingredients, such as the stellar atmosphere calcula¬ 
tions, are nevertheless affected by their own uncertainties). 


Fig-Dshows that, regardless of the methodology, the Pb 
exhibits a far more dispersed correlation with [Mg/Fe] than 
no, i.e. [Mg/Fe] is not the main driver of the IMF variations 
in ETGs. The fact that Tb tends to marginally increase with 
[Mg/Fe] mostly stems from the lack of galaxies with low- 
[Mg/Fe] (< 0.1) at the highest ao (see red symbols in the 
right panels of Fig. |T]) and the fact that [Mg/Fe] increases 
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Table 1. Methods to derive the IMF slope, F^,. Col. 1 labels the method, while cols. 2 and 3 report the selected list of fitted spectral 
features and elemental abundances ([X/Fe]), respectively. For the lSSP+[X/Fe] approach, we include the 1SSP features (from LB13), plus 
other Lick-based abundance-sensitive indices (marked with superscripts in col. 2). For each [X/Fe] in col. 3, we report the superscripts 
corresponding to the spectral feature(s) with prominent sensitivity to it. 


Method 

(1) 

EW 

(2) 

[X/Fe] 

(3) 

1SSP 

H/3 0 , [MgFe]', TiOl, Ti02, Mg4780, Na8190 SD ss, Ca2, CaHK, NaD 

none 

lSSP+[X/Fe] 

H/3 0 , [MgFe]', TiOl, Ti02, Mg4780, Na8190 SD ss, Ca2, Cal, Ca4227 x , 
NaD 2 , Fe4531 3 , Mgl 4 , Mg2 5 , C4668 6 , CN2 7 , Mgb5177 8 

[Ca/Fe] 1 , [Na/Fe] 2 , [Ti/Fe] 3 , 
[O/Fe] 4 ’ 7 , [C/Fe] 6 ’ 4 
[N/Fe] 7 , [Mg/Fe] 5 - 8 , [Si/Fe] 4 ’ 5 


Table 2. Best-fitting coefficients, Ai, A2, and A3 (cols. 2-4), of 
Eq.H for different methods to derive T± (col. 1). The row labeled 
“• • -±[Mg/Fe]*” corresponds to results for lSSP+[X/Fe] fits when 
replacing [Mg/Fe] in Eq.[2] with [Mg/Fe]* (see the text). 


Method 

(1) 

Ar 

(2) 

A2 

(3) 

A 3 

(4) 

1SSP 

3.76 ±0.16 

-0.08 ±0.40 

1.95 ±0.02 

lSSP+[X/Fe] 

2.92 ±0.24 

-0.93 ±0.15 

1.60 ±0.02 

■ ■ ■ + [Mg/Fe]* 

2.80 ±0.22 

-0.90 ±0.14 

1.56 ±0.02 


with <jo in ETGs (e.g. lThomas et ai1l2005l ). To illustrate this 
point, we fit a bivariate relation to the data, 

T b = Ai logcro+A 2 [Mg/Fe] + A 3 . (2) 

This approach is the same as ISmi ths (i2014l . hereafter S14), 
with the noticeable difference that in the present work we 
study the IMF slope, Tb, rather than its overall normal¬ 
ization. Eq. [2] allows us to single out, quantitatively, the 
contribution of [Mg/Fe] and a o to the systematic varia¬ 
tions in the IMF. We derive Ai, A 2 , and A 3 (see Tab. EJ 
with a robust least-square fitting procedure, minimizing ab¬ 
solute residuals to Tb, excluding those stacks with large 
(> 50%) uncertainty on T b (see grey symbols in Fig. |T]). 
Fig. El plots the best-fitting slopes, Ai and A 2 , for different 
fitting methods together with their bootstrapped confidence 
contours, with red and blue colours corresponding to 1SSP 
and lSSP+[X/Fe] results, respectively. The green dot and 
ellipses represent the lSSP+[X/Fe] slopes obtained by re¬ 
placing [Mg/Fe] in Eq. El with [Mg/Fe]*, i.e. the best-fitting 
[Mg/Fe] from lSSP+[X/Fe] fits. In the 1SSP method, no 
significant correlation with [Mg/Fe] is found at fixed 00 , as 
A 2 is consistent with zero at less than 1 a, whereas in the 
lSSP±[X/Fe] case, A 2 is significantly negative, i.e. F b tends 
to decrease with [Mg/Fe] at fixed 00 . Notice that between 
these two methods, there is a ~ 20% difference in Ai, with 
shallower slopes for the lSSP+[X/Fe] case, consistent with 
Fig.m Furthermore, given the range of logon (~ 1.2 dex) 
and [Mg/Fe] (~ 0.3 dex) probed by the data, the values of 
Ai and A 2 for the lSSP+[X/Fe] method imply a narrower 
variation in the IMF slope at fixed velocity dispersion, with 
respect to the variations associated to 00 . Hence, both ap¬ 
proaches agree on [Mg/Fe] having a minor role in the IMF 
trend of ETGs. 


5 SUMMARY AND DISCUSSION 

We have assessed in this letter the relative role of central 
velocity dispersion and [Mg/Fe] abundance ratio as drivers 
of the systematic variation found in the IMF of ETGs. We 
assemble a set of 126 high-quality spectra from a sample of 
24,781 nearby ETGs (LB13). We split the range in central 
velocity dispersion (100<ao<320 kms -1 ) into 18 bins. For 
each bin, i.e. at fixed 00 , we create six stacks according to 
[Mg/Fe], covering a range A[Mg/Fe]~ 0.3dex. For each of 
the 126 stacks (i.e. including the original 18 stacks in LB13 
only binned with respect to on), we derive the slope of the 
(bimodal) IMF, Tb, assuming a single burst, and considering 
both a direct fit to the line strengths (1SSP method, as in 
LB13), and an alternative approach (lSSP+[X/Fe[), where 
IMF and individual abundance ratios are fitted simultane¬ 
ously. 

- We find a very mild variation of Fb with [Mg/Fe] at fixed 
(Jo, implying that [Mg/Fe] plays a subdominant role in the 
variation of the IMF slope among ETGs. 

- On the contrary, we find a strong, highly significant (> 
10 ct), correlation of the IMF slope with 00 , regardless of the 
adopted methodology. 

Since for a fixed IMF functional form (such as the bimodal 
distribution), Tb is proportional to the overall normalization 
(i.e. the mass-to-light ratio normalized to that for a Kroupa- 
like distribution), our results contrast with the claim of 
CvD12b, that the IMF should correlate more strongly with 
[Mg/Fe] than oq. We note that the models of CvD12b use 
[Fe/H] and not [Z/H] as an “overall” metallicity parameter. 
Therefore their [Mg/Fe] variations could be related to to¬ 
tal metallicity. S14 found that for a sample of 34 ETGs in 
common between CvD12b and Atlas 3D , the dynamical con¬ 
straints to the IMF appeared not to be consistent with the 
line strength analysis of CvD12b. Dynamical models cor¬ 
relate with velocity dispersion (cro) and not with [Mg/Fe]. 
Therefore, our results - based on spectral indices - are qual¬ 
itatively more consistent with the results from dynamical 
modeling. In particular, the fact that Tb decreases slightly 
with [Mg/Fe], at fixed cro, when using the lSSP±[X/Fe] fit¬ 
ting method, is similar to the results in S14 regarding the 
dynamical constraints (see fig. 4 of S14, and our Fig. El- 
At the moment, the origin of the above discrepancies is 
not clear. As already noticed by S14, several possible effects 
might contribute, including (i) sample size/selection issues, 
as we analyze a sample of 126 stacked spectra, while CvD12b 
studied 38 (individual) galaxies; (ii) methodology, as we rely 


© 2015 RAS, MNRAS 000. [T][6] 














No correlation of IMF with [Mg/Fe] 5 



Figure 2. Best-fitting slopes of the bivariate correlation of Tb 
vs. log cro and [Mg/Fe] (Eq. (2j- Different methods/samples are 
plotted with different colors, as labeled in the upper-right of the 
figure. Ellipses mark 1 cr and 2 a normal confidence contours on 
the slope values. The horizontal dashed line represents a zero cor¬ 
relation with [Mg/Fe] at fixed cro. Notice that the log cro slope is 
highly significant in all cases, while the [Mg/Fe] slope depends on 
the adopted methodology, being fully consistent with zero when 
no abundance ratios are fitted to spectral indices (red dot and 
ellipses), or mildly anticorrelated if individual abundances are 
considered. 


on line strengths, rather than spectral fitting, and adopt 
(MIUSCAT) stellar population models with varying total 
metallicity (in contrast to CvD12a models); (iii) systematics, 
on either stellar population models (for spectroscopic con¬ 
straints) or the modeling of the (dark-)matter distribution 
in ETGs (for dynamical constraints); (iv) radial variations 
of the IMF, as the SDSS spectra cover ~ lR e ff spatial scale, 
in contrast to the R e ff/8 scale of CvD12b. The latter aspect 
is particularly relevant in light of our recent finding (MN15) 
that the IMF slope decrea ses with galacto-centr ic distance 
in massive ETGs (see also IPastorello et al.ll2014r ). implying 
that local and global constraints to the IMF might result in 
different trends. We also remark that the present work does 
not imply that the central no (an integrated measurement) 
is the “causal” driver of the IMF in ETGs, but, rather, that 
[Mg/Fe] (possibly related to the timescale of star formation) 
is not driving the systematic variations of the IMF, at least 
over a ~ lR e ff spatial scale. In particular, given the tight 
correlation between metallicity and oo in ETGs, and the 
narrow dynamical range of metallicity in our stacked spec¬ 
tra at fixed o o (LB13) , the values of Tb in Fig. [l]can also 
be considered to correlate with (total) metallicity. Hence, 
our integrated spectra cannot be used to infer the physical 
driver of the trend. Spatially resolved spectroscopic studies, 
extended to large samples of galaxies, with the aid of further 
advances in the modeling of stellar populations, should help 
in the future to address this issue (Martfn-Navarro et al., in 
preparation). 


ACKNOWLEDGMENTS 

We thank Dr. M. Beasley, and the anonymous 
referee, for interesting comments and sugges¬ 
tions. We acknowledge the use of SDSS data 
(http://www.sdss.org/collaboration/credits.html), and 
support from grant AYA2013-48226-C3-1-P from the Span¬ 
ish Ministry of Economy and Competitiveness (MINECO). 


REFERENCES 

Adelman-McCarthy, J. K., Agiieros, M. A., Allam, S. S., 
et ah, 2008, ApJS, 175, 297 
Auger, M. W., et ah, 2010, ApJ, 724, 511 
Cappellari, M., et ah, 2012, Nature, 484, 485 
Cappellari, M., et ah, 2013, MNRAS, 432, 1862 
Chabrier, G., PASP, 115, 763 

Chabrier, G., Hennebelle, P., Chariot, S., 2014, ApJ, 796, 
75 

Cid Fernandes, R., Mateus, A., Sodre, L., Stasinska, G., 
Gomes, J. M., 2005, MNRAS, 358, 363 
Cenarro, A. J., Gorgas, J., Vazdekis, A., Cardiel, N., 
Peletier, R. F., 2003, MNRAS, 339, L12 
Conroy, C., van Dokkum, P., 2012a, ApJ, 747, 69 (CvD12a) 
Conroy, C., van Dokkum, P., 2012b, ApJ, 760, 71 (CvD12b) 
de La Rosa, I. G., La Barbera, F., Ferreras, I., de Carvalho, 
R. R., 2011, MNRAS, 418, L74 
Dutton, A. A., Maccio, A. V., Mendel, J. T., Simard, L., 

2013, MNRAS, 432, 2496 

Ferreras, I., La Barbera, F., de la Rosa, I. G., Vazdekis, 
A., de Carvalho, R. R., Falcon-Barroso, J., Ricciardclli, 
E„ 2013, MNRAS, 429, L15 (F13) 

Ferreras, I., Weidner, C., Vazdekis, A., La Barbera, F., 
2015, MNRAS Lett, in press (arXiv:1501.01636) 

Hopkins, P. F., 2013, MNRAS, 433, 170 
Johansson, J., Thomas, D., Maraston, C., 2012, MNRAS, 
421, 1908 

Kroupa P., Weidner C., Pflamm-Altenburg J., et ah, 2013, 
in Vol. 5: Galactic Structure and Stellar Populations, ed. 
T. D. Oswalt & G. Gilmore (New York: Springer), 115 
Kuntschner, H., et ah, 2010, MNRAS, 408, 97 
La Barbera, F., de Carvalho, R. R., de la Rosa, I. G., Lopes, 
P. A. A., Kohl-Moreira, J. L., Capelato, H. V., 2010a, 
MNRAS, 408, 1313 (Paper I) 

La Barbera, F., et ah, 2013, MNRAS, 433, 3017 
McDermid, R. M., et ah, 2014, ApJ, 792, 37 
Martfn-Navarro, I., La Barbera, F., Vazdekis, A., Falcon- 
Barroso, J., Ferreras, I., 2015, MNRAS, 447, 1037 (MN15) 
Pastorello, V., et ah, 2014, MNRAS, 442, 1003 
Sansom, A. E., de Castro Milone, A., Vazdekis, A., 
Sanchez-Blazquez, P., 2013, MNRAS, 435, 952 
Smith, R., Lucey, J. R., 2013, MNRAS, 434, 1964 
Smith, R., MNRAS, 443, 69 (S14) 

Spiniello, C., Trager, S., Koopmans, L. V. E., Conroy, C., 

2014, MNRAS, 438, 1483 (STK14) 

Spiniello, C., Trager, S., Koopmans, L. V. E., 2015, 
(arXiv: 1401.1145) 

Thomas, D., et ah 2005, ApJ, 621, 673 
Thomas, D., Maraston, C., Johansson, J., 2011, MNRAS, 
412, 2183 

Tortora, C., Romanowsky, A. J., Napolitano, N. R., 2013, 
ApJ, 765, 8 


© 2015 RAS, MNRAS 000, [GEO 








6 F. La Barbera, I. Ferreras, A. Vazdekis 


van Dokkum, P. G., Conroy, C., 2010, Nature, 468, 940 
Vazdekis, A., Casuso, E., Peletier, R. F., Beckman, J. E., 
1996, ApJS, 106, 307 

Vazdekis, A., Ricciardelli, E., Cenarro, A. J., Rivero- 
Gonzalez, J. G., Diaz-Garca, L. A., Falcon-Barroso, J., 
2012, MNRAS, 424, 157 (MIUSCAT-I) 


© 2015 RAS, MNRAS 000, [Hi] 


